AD-A134  Oil  A  POLARIZATION  DIVERSITY  RADAR  DATA  PROCESSOR ( U )  AIR 

FORCE  GEOPHYSICS  LAB  HANSCOM  AFB  MA  d  I  METCALF  ET  AL 
28  APR  83  AFGL  TR-83-01 1  I 


UNCLASSIFIED 


F/G  17/9 


I 


► 


MICROCOPY  RESOLUTION  TEST  CHART 
NinoMi  mui  V 


'  'v  •  .*•  -.4  '*  •  .  “  -•  . 


Unclassified 


SECURITY  CLASSIFICATION  QW  THIS  PAGE  (Whan  ggAggfAFAj) 


REPORT  DOCUMENTATION  PAGE 

READ  mrrRUCTtOMS 

BEFORE  COMPLETING  FORM 

1.  HC.O.T  NUM.E*  ].  OOVT  ACCESSION  NO. 

AFGL-TR-83-0111 

1-  RECIPIENT'S  CATALOG  NUMBER 

ll 

4.  TITLE  !md  Subltll.)  ' 

A  POLARIZATION  DIVERSITY  RADAR 

DATA  PROCESSOR 

5  TYPE  OF  REPORT  A  PERIOD  COVERED 

Scientific.  Interim. 

S-  PERFORMING  OTG.  REPORT  NUMRER 

IP,  No.  318 

7.  AUTHONfa; 

James  I.  Metcalf 

Graham  M.  Armstrong 

B.  CONTRACT  OR  GRANT  NUMBER^; 

t.  performing  organization  name  and  address 

Air  Force  Geophysics  Laboratory  (LYR) 

Hans  com  AFB 

Massachusetts  01731 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMRERS 

62101F 

66701602 

II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Air  Force  Geophysics  Laboratory  (LYR1 

Hanscom  AFB 

Massachusetts  01731 

t2.  REPORT  DATE 

28  April  1983 

1!  NUMRER  OP  PAGES 

27 

M.  MONITORING  AGENCY  NAME  A  ADDRESS (II  dlllarmnt  from  Controtllni  Olftca) 

IS  SECURITY  CLASS,  (of  thla  rapori) 

Unclassified 

IS*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

»«.  DISTRIBUTION  STATEMENT  (ol  thla  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  tha  abstract  antarad  in  Block  20.  If  different  front  Rap ort) 


supplementary  notes 


19.  KEY  WORDS  (Conllnu*  on  rovmrmo  mlda  If  nacaaaary  And  Identify  by  block  nuoib*r) 

Polarization-diversity  weather  radar 
Radar  data  processor 
Polarization  switching 
Signal  autocovariance 


40  ABSTRACT  ( Conllnua  on  ravaraa  tldm  II  nfeiMAT  and  Identify  by  block  nmmbar) 

^  A  real  time  data  processor  has  been  designed  for  use  with  the  AFGL 
10- cm  Doppler  weather  radar  which  is  to  be  operated  with  alternating  trans¬ 
mission  of  horizontally  and  vertically  polarized  signals.  In  this  mode  of  opera¬ 
tion  the  reception  of  backscattered  signals  of  polarizations  identical  to  those  of 
the  transmitted  signals  allows  the  computation  of  the  differential  reflectivity 
between  the  two  polarizations  in  addition  to  the  absolute  reflectivity  and  the 
Doppler  mean  velocity  and  variance.  The  switching  of  transmitted  .polarization 
_introduces^difficultiesjjT_the_estimation_of>Jhe_autocovariance_of>the>£eceilved^; 

00  I  JAN  7}  1473  coition  or  i  nov  ts  is  obsolete  Unclassified 


SECURITY  CLASSIFICATION  OP  THIS  PAGctfPhAft  Bnlarod) 


Unclassified 


20.  (Contd) 

'signals,  from  which  the  Doppler  velocity  parameters  are  derived.  The 
processor  design  and  the  allowed  modes  of  radar  operation  circumvent  these 
difficulties.  This  report  describes  the  processing  algorithms  theoretically 
and  presents  details  of  the  implementation  of  these  algorithms  in  hardware. 


Accession  For 

NTI3  GRA&I  W 

DTIC  TAB 

□ 

Unannounced 
Justification _ 

□ 

[  Bv _ 

^Distribution/  _ 

Availability  Codes 
Avail  and/or 
Dlst  Special 


m 


Unclassified 

JtCUSITY  CL  *»»iriC»TIO«  D««  KnlfMi 


IF 


Contents 


1.  INTRODUCTION  5 

2.  SIGNAL  PROCESSING  REQUIREMENTS  6 

2.  1  Doppler  Velocity  Parameters  6 

2.  2  Absolute  and  Differential  Reflectivity  13 

3.  SIGNAL  PROCESSOR  DESCRIPTION  15 

3.  1  Autocovariance  Calculator  17 

3.2  Average  Velocity  Calculator  19 

3.  3  Doppler  Variance  Calculator  20 

3.4  Reflectivity  Calculator  21 

3.  5  Timing  Generator  23 

4.  OPERATIONAL  PROCEDURES  AND  PARAMETERS  23 

5.  CONCLUSIONS  26 

REFERENCES  27 


Illustrations 


1.  Standard  Deviation  of  Estimate  of  Differential  Reflectivity  15 

2.  Block  Diagram  of  the  Radar  Data  Processor  16 

3.  Details  of  the  Autocovariance  Calculator  18 

4.  Diagram  of  Mean  Velocity  Calculator  19 


Illustrations 


5.  Details  of  the  Doppler  Variance  Calculator  20 

6.  The  Reflectivity  and  Differential  Reflectivity  Calculator  22 

7.  Block  Diagram  of  Timing  Generator  24 

8.  Timing  Diagram  25 


Tables 


1.  Value  of  Correction  Index  for  Use  With  Autocovariance  at 

Two-pulse  Lag  7 

2.  Percentage  of  Pulse  Pairs  Usable  for  Autocovariance  Calculation  13 

3.  Data  Processor  Characteristics  25 


4 


A  Polarization  Diversity  Radar  Data  Processor 


1.  INTRODUCTION 

The  development  of  a  polarization  diversity  capability  for  the  AFGL  10-cm 
1  2 

Doppler  weather  radar  '  necessitates  the  design  and  construction  of  a  new  data 
processor  to  compute  reflectivity,  Doppler  parameters,  and  polarization  parameters 
in  real  time.  The  first  phase  of  radar  modification  required  for  implementation  of 
this  capability  was  begun  in  1982,  with  the  long  range  goals  of  transmitting  either 
circularly  or  linearly  polarized  signals,  with  rapid  switching  between  either  right 
and  left  circular  or  horizontal  and  vertical  polarizations,  and  coherently  receiving 
two  signals  of  polarizations  identical  and  orthogonal  to  that  of  the  transmitted  signal. 
The  first  phase  involves  modification  of  the  antenna  and  feed  system  to  a  dual  polari¬ 
zation  configuration  and  the  installation  of  a  high  power  microwave  switch.  When 
this  phase  is  completed,  as  is  planned  for  FY  1984,  the  radar  will  have  the  capability 
of  transmitting  horizontally  and  vertically  polarized  signals,  with  rapid  switching 
between  the  two  polarizations,  and  receiving  backscattered  signals  of  polarization 


(Received  for  publication  26  April  1983> 
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identical  to  that  of  the  transmitted  signals.  The  radar  can  then  be  used  to  measure 

3 

the  differential  reflectivity  as  described  by  Seliga  and  Bringi.  This  quantity, 
together  with  the  absolute  reflectivity  derived  from  backscattered  signals  of  one 
polarization,  provides  information  on  the  shapes  of  particles  in  the  backscattering 
medium  which  can  be  interpreted  in  terms  of  thermodynamic  phase  and,  in  rain, 
drop  size.  The  measurement  and  interpretation  of  differential  reflectivity  has  been 
discussed  extensively  in  the  scientific  literature. 

The  Doppler  mean  velocity  and  variance  derived  from  coherently  received 
backscatter  signals  yield  further  information  on  the  backscattering  medium  and 
provide  a  measure  of  the  air  flow  patterns  in  precipitation  systems.  Taken  to¬ 
gether  with  the  absolute  and  differential  reflectivities,  the  Doppler  velocity  param¬ 
eters  permit  the  derivation  of  a  comprehensive  and  detailed  description  of  the  pre¬ 
cipitation  system. 

The  concepts  and  formulations  of  the  data  processing  are  described  in  the 
following  section.  Section  3  describes  the  implementation  of  the  algorithms  and 
includes  block  diagrams  of  the  processor  and  its  components. 

2.  SIGN  AL  PROCESSING  REQUIREMENTS 

In  this  section  we  describe  the  requirements  and  procedures  for  processing 
Doppler  velocity  data,  derived  from  the  coherent  receiver,  and  reflectivity  data, 
derived  from  the  logarithmic  receiver.  Procedures  are  described  for  the  two  cases 
of  radar  operation  with  and  without  polarization  switching.  Of  particular  concern 
for  both  reflectivity  and  Doppler  velocity  computations  are  the  effects  of  polariza¬ 
tion  switching  on  the  accuracy  and  possible  bias  of  the  computed  quantities.  Hence, 
we  discuss  several  problems  inherent  in  operation  with  switched  polarization  and 
describe  the  wavs  in  which  these  problems  are  solved  or  circumvented  in  the  pro¬ 
cessor  design  and  radar  operation. 

2. 1  Doppler  Velocity  Parameter* 

The  Doppler  mean  velocity  and  variance  are  to  be  computed  by  means  of  pulse 
pair  autocovariance  techniques.  In  particular,  the  mean  velocity  v  is  to  be  calculated 

4 

from  the  autocovariance  at  one  and  two  pulse  lage  using  the  technique  due  to  Rummler 


3,  Seliga,  T.  A.  ,  and  Bringi,  V.  N.  (1976)  Potential  use  of  radar  differential 

reflectivity  measurements  at  orthogonal  polarizations  for  measuring 
precipitation,  J.  Appl.  Meteor.  ,1^5 :^6 9-76. 

4.  Rummler,  W.D.  (1968)  Introduction  of  a  New  Estimator  for  Velocity  Spectral 

Parameters,  Tech.  Memo.  MR1-68-4  121-5,  Sell  Telephone  Laboratories, 
Whippany,  New  Jersey. 
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(1) 


=  1  '  Wt  arg  lR(T  *1 


and  the  technique  suggested  by  Strauch  et  al. 
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v  ”  7  '  TiTf  arg  [R(2t  )]  +  n7T] 


(2) 


where  the  value  of  the  correction  index  n  is  determined  by  the  arguments  of  R(T  1 
and  R(2t)  as  shown  in  Table  1. 


Table  1.  Value  of  Correction  Index  for  Use  With 
Autocovariance  at  Tivo-pulse  Lag.  Value  is 
determined  by  the  quadrants  of  the  arguments  of 
R(t  )  and  R(2r  1.  Brackets  denote  conditions  for 
which  derived  velocity  is  negative 


arg  [  R(T  1] 

(-17,  -77/2) 

arg  [R(2r  1] 

(-77/2.  01  (0.  77/2) 

(il/ 2,  77) 

CVJ 
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1 

1 

[0] 

1 

1-1] 

1-1] 

(-7T/2,  0) 

[0] 

(0] 

0 

1-1] 

(0.  77/2) 

1 

10] 

0 

0 

(it  12,  n) 

1 

1 

1-1) 

0 

Several  approaches  exist  for  the  estimation  of  Doppler  velocity  variance.  If 
the  spectrum  is  Gaussian  and  the  noise  level  is  known,  then  an  exact  formula  for 
the  variance  is 


TTT 

877  T 


In 


R(o)  -  N 

~  "RTrT" 


(3) 


where  N  is  the  average  noise  level. 


The  formula  developed  by  Rummler 
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is 
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I  R(T  )  I 
R(o) 


(1  +  SNR 


’’I 


(4) 


5.  Strauch,  R.  G. ,  Kropfli,  R.  A.,  Sweezy,  W.  B. ,  Moninger,  W.  R. ,  and 

Lee,  R.  W.  (1978)  Improved  Doppler  velocity  estimates  by  the  poly-pulse- 
pair  method.  Preprints,  18th  Conf.  Radar  Meteor. ,  Amer.  Meteor.  Soc. , 
Boston:  pp.  376-380. 
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where  SNR  is  the  signal-to-noise  ratio.  The  noise  correction  can  be  omitted  in 
either  of  the  above  formulae  if  the  signal-to-noise  ratio  is  high.  A  preferred 

g 

alternative  is  the  formula  proposed  by  Srivastava  et  al. 


24/7  2t  2 


In 


R(T  ) 
T?T2t7 


which  does  not  require  real-time  estimation  of  the  signal-to-noise  ratio.  The 
formula  is  also  exact  for  Gaussian  spectra.  Other  estimation  procedures  yield 
better  results  if  the  Doppler  spectrum  is  significantly  non-Gaussian.  '  These  are 
all  based  on  the  computation  of  the  autocovariance  of  the  received  signal  at  several 
time  lags.  Hence,  although  we  do  not  contemplate  implementing  these  procedures, 
we  note  that  they  can  be  easily  added  as  options  in  the  future,  since  the  present 
processor  design  includes  the  real-time  computation  of  autocovariance  at  zero, 
one,  and  two-pulse  lags. 

The  key  problem  in  the  processor  design  is  therefore  to  implement  optimum 
estimators  of  the  autocovariance  function,  that  is,  estimators  that  are  unbiased  in 
both  magnitude  and  phase.  The  autocovariance  at  zero  time  lag  will  be  a  sum  of 
two  categories  of  terms,  corresponding  to  horizontal  and  vertical  transmitted 
polarization,  respectively.  For  complex  amplitudes  AH(0nj  and  Av(2n+  at 
horizontal  (HI  and  vertical  (VI  polarizations  and  successive  times  designated  by 
2n  and  2n  +  1,  we  have 

{-1 

K<0>  =  w  {  E  AH(2nl)AH 
v  n  =  u 

=  1/2  |  RHh(0>+  *W0> 


2n) 


N  i 

T  1 


n=  0 


AV(2n+  l)AV(2n 


+  1>i 


(6) 


where  Rjj^(O)  and  Ryy(O)  are  the  estimates  one  would  derive  separately  from  sig¬ 
nals  of  horizontal  and  vertical  polarization,  respectively.  The  expected  values  of 
these  estimators  are  the  average  received  powers  for  the  two  polarizations,  namely 


6.  Srivastava,  R.  C.  .  Jameson,  A.  R. ,  and  Hildebrand,  P,  H.  (19791  Time-domain 

computation  of  mean  and  variance  of  Doppler  spectra.  J.  Appl.  Meteor, 
jyj,:  189-194.  - 

7.  Passarelli,  R.  E. ,  Jr.,  and  Siggia,  A.  D.  ( 1 983)  The  autocorrelation  function 

and  Doppler  spectral  moments:  Geometric  and  asymptotic  interpretations. 

J.  Clim.  Appl.  Meteor.  2g_,  accepted  for  publication 
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(7) 
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e[Rhh(0>]  = 

E|Rvv(0>]  = 

For  further  analysis  of  the  autocovariance  estimators  one  can  normalize  the  signal 

1/2  1/2 

amplitudes  by  their  root-mean-square  values,  (Fpj)  '  and  (F,.)  .  The  bias  of 

K(0)  relative  to  the  estimate  obtained  from  a  radar  operating  with  fixed  polarization, 
for  example,  horizontal,  is  determined  by  the  ratio  of  received  power  between  the 
two  polarizations,  that  is,  the  differential  reflectivity.  Thus-. 

fc[HHH(0)J  +  K(Kvv(0)] 

=  d/2)  |eihhh«»]  ElrWo)]J 

=  (1/2)  (  1  +Fv/F(1)E|Rhh(0)]  .  (8) 


E[H(0)J  =  (1/2) 


Computation  of  the  autocovariance  at  non-zero  time  lags  requires  evaluation  of 
four  categories  of  signal  pulse  pairs:  (1)  those  having  common  horizontal  polariza¬ 
tion  of  the  two  transmitted  pulses,  (2)  those  having  common  vertical  polarization, 
(3)  those  in  which  the  first  pulse  is  horizontally  polarized  and  the  second  vertically 
polarized,  and  (4)  those  in  which  the  first  pulse  is  vertically  polarized  and  the 
second  horizontally  polarized.  We  can  then  analyze  the  most  general  case  in  which 
any  number  of  pulses  may  be  transmitted  during  the  polarization  switching  interval. 
For  analysis  of  the  autocovariance,  we  represent  phase  of  the  received  signals 
explicitly,  in  terms  of  phase  due  to  scattering  and  phase  due  to  signal  propagation 
through  the  radar  and  the  atmosphere.  Because  the  transmitted  signals  of  ortho¬ 
gonal  polarization  follow  different  paths  in  the  radar  and  because  precipitation  is 
typically  an  anisotropic  propagation  medium,  there  is  a  nearly  constant  differential 
phase  shift  (in  a  given  measurement  scenario)  between  received  signals  corres¬ 
ponding  to  orthogonally  polarized  transmissions.  Thus  we  have  for  horizontal 
polarization 


II 


=  A 


i(0g<  T)  + 


V 


H1 


(9) 


and  for  vertical  polarization 

,  ,  i<os(  T)  +  V 

Av  =  ^Av^  e 


(10) 


J 

i. 
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where  9s(t)  is  the  phase  due  to  scattering  (assumed  to  be  independent  of  polariza¬ 
tion)  and  and  9y  are  phases  due  to  propagation  in  the  radar  and  in  the  atmosphere 
at  horizontal  and  vertical  polarization,  respectively.  The  four  categories  of  terms 
available  for  the  autocovariance  estimator  R(7  )  are  represented  by 


|  Ah(0)Ah(t  )  = 

|ah(0)| 

1  Ah(t  )| 

i[®  (0)  -  ©  (r)] 
e  s  s 

(11) 

|  Av(0)Av(r  )  = 

! 

Iav(0)| 

Iav(t  )l 

i{<p  (0)  -  ©  (t)] 
e  s  s 

(12) 

Ah(0)A*(t)  = 

i 

1 A  H(  0)  1 

1  Ay(r  )| 

i[((?s(0)  -  9s(t))+  (9h  -  9V)] 

(13) 

I 

1  A,.(0)Au(T)  = 

1 

I  A v(0>  | 

!ah(t  >| 

i[(9g(0)  -  9s(r))  -  (9h  -  9V)J 

(14) 

We  assume  that  the  phase  difference  0  (0)  -  0  (r  )  =  A<?(t  )  is  a  random  variable 

=  s 


independent  of  polarization  and  independent  of  time.  The  sums  of  the  four  cate¬ 


gories  of  pulse  pairs  correspond  to  autocovariance  estimators  RHjj(7  ),  R^.(7  ), 


R^.(T  ),  and  Rv(j(T  ).  By  scaling  the  amplitudes  as  described  above  we  can  ex¬ 


press  the  expected  values  of  three  of  these  estimators  in  terms  of  the  expected 


value  of  Kjjj.j(T  ): 


E(RVV(T  )]  =  <FV/PU>  E|Rhh(t  )] 

,,,  -  i(0  -  0  ) 

E|R(iv(T  )]  =  fPv/r>H>1'^  E|Rmi(T  >]  e 


(15) 


(16) 


-  1')  "  i(J  1 1  "  d  .  .  ) 

E[RVIi(T>]  =  <Pv/Fn>  “  E[R|ui(t  )]  e 


(17) 


An  analogous  derivation  can  be  developed  for  the  computation  of  R(2r  ).  From  the 
results  above,  we  can  develop  three  estimators  of  R(T  ).  (  sing  equal  numbers  of 

pairs  of  common  horizontal  and  vertical  polarization  we  have 


Rj(T)  =  (1  ' 2)  I  RIItI(T  )  -  KVv(T‘]  • 


(18) 


The  expected  value  of  this  estimator 

E[K,(7)]=  ( 1/2)  ( 1  +Pv/F|f]  E(  R[1h<t  >] 


(19) 
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is  unbiased  in  phase  but  is  biased  in  magnitude  by  a  function  of  the  differential 
reflectivity.  A  second  estimator  involves  pairs  of  the  dissimilar  polarizations, 
in  equal  numbers  of  those  with  horizontal  first  and  those  with  vertical  first: 


R2(t)  =  (1/2)|Rhg(t)  +  Rvh<t>]  . 


(20) 


The  expected  value  of  this  estimator 

E|R2(t)]  =  (Trv/Fa)1/2  E|Rhh(t)]  cosW>H  -  0V) 


(21) 


is  again  unbiased  in  phase  but  is  biased  in  magnitude  by  the  square  root  of  the  power 
ratio  and  by  the  cosine  of  the  differential  phase  shift.  The  differential  phase  shift 
in  the  system  can  be  minimized  by  careful  design.  However,  the  differential  phase 
shift  in  the  atmosphere  increases  monotonieally  with  range,  and  at  long  range 
through  heavy  precipitation  it  can  approach  90°,  at  which  point  the  expected  value 
of  R^tT  )  is  zero.  This  estimator  could  be  used  for  mean  velocity  estimation,  by 
Eq.  (1),  but  would  yield  highly  biased  variance  estimates.  A  third  estimator  can  be 
developed  by  recognizing  that 


E|KI1V<T  )] 

t:n\  „<t  >1 


1/2 


i(d. 


V 


(22) 


■'*  *  1  /  •) 

We  use  the  ratio  |  R^.IT  *1  “as  an  estimator  of  the  phase  factor 

exp|i(0(J  -  0^  )]  and  define  the  estimator  as 


R.j(T  ) 


(«hv(t» 

■Rvn<r)' 

R,„.<T  ) 

\ 

L  uv 

1/2 


+  KV1I(T) 


»H\  (T) 

«vh<7 ) 


1/2 


khv(t  '  rvh(t  * 


1/2 


(23) 


This  estimator  would  be  unbiased  only  if  its  expected  value  were  equal  to  the  product 
of  the  expected  values  of  |  Rjj^.(T  )]'  2  and  1  R^.jj(T  )]' ^2,  that  is,  if  Rj^.lT  )  and 

jj(r  )  were  uncorrelated.  Typically,  however,  these  estimates  are  highly  corre¬ 
lated,  especially  if  they  are  derived  from  the  same  sample  time  series.  In  general, 
for  correlated  random  variables  x  and  v  we  have 


E(xv)  =  E(x)E(v)  +  Cov  (x,y) 


(24) 


and 


Cov  (x,  y)  =  p(x,  y)  [Var(x)  Varlvl] 


where  Cov  is  the  covariance,  p  the  autocorrelation,  and  Var  the  variance.  In  the 
present  case,  p  ;  1  so  that  if  x-  !  (T  )] 1  ^  and  y  =  |  R^. j^lr  )] 1 '  2 


E[R3(t)]  -  E[[Rhv(t  )]1/21  E|[Rvu(t  ))1/2] 

+  [Var[(RHV(T))1/2]  Var|afVH(T))1/2])1/2  . 


VV'e  have  not  fully  analyzed  the  statistical  properties  of  this  estimator.  The  vari¬ 
ances  are  approximately  inversely  proportional  to  the  number  of  pulse  pairs  in  the 
estimators  R...,(T  )  and  R,.0(T  ),  and  can  therefore  be  reduced  bv  extending  the 
length  of  the  signal  sample.  In  practice,  the  available  length  of  the  sample  is 
limited  by  constraints  on  scanning  time  and  spatial  resolution.  Hence,  it  is  not 
certain  that  the  variances  can  be  made  negligibly  small  or  even  "acceptably"  small. 

In  summary,  we  find  that  reliable  estimation  of  autocovariance  can  be  accom¬ 
plished  only  by  averaging  the  products  of  signals  having  common  polarization.  The 
bias  factor  ( 1  +  which  arises  from  the  collective  summation  of  pairs  of 

horizontally  polarized  signals  and  pairs  of  vertically  polarized  signals  does  not 
affect  the  ratio  of  autocovariances  at  different  time  lags,  as  used  in  Eq.  (5). 

The  problems  inherent  in  the  use  of  summations  of  products  of  signals  of 
opposite  polarizations  can  be  alleviated  in  our  system,  as  the  Doppler  processing 
and  reflectivity  processing  are  done  with  signals  transmitted  on  different  carrier 
frequencies  which  can  be  pulsed  at  different  rates.  Specifically,  the  carrier  fre¬ 
quency  for  the  coherent  receiver  can  be  pulsed  at  1,  2,  3,  or  4  times  the  pulse 
repetition  rate  of  the  carrier  frequency  for  the  logarithmic  receiver.  This  feature 
was  included  in  the  original  radar  design  for  the  purpose  of  resolving  range  and 
velocity  ambiguities.  If  the  multiples  3  or  4  are  used,  then  sums  of  products  of 
signals  of  common  polarization  can  be  accumulated  to  estimate  the  autocovariance 
at  1  and  2  pulse  time  lags.  Thus  we  will  implement  the  estimators 

R1(r)=  (1/2>(Rhh(t)  +  Rvv<r)J  (27) 


R1(2t)  =  (1/2)[Rhh(2t)  +  Rvv(2r)]. 


•;*r,  -‘t 


The  switching  of  transmitted  polarization  will  signal  the  processor  to  delete  from 
the  summation  those  pairs  of  pulses  spanning  the  change  of  polarization.  Table  2 
shows  the  percentage  of  all  pairs  that  are  usable  for  the  autocovariance  calculation, 
that  is,  that  have  common  polarization  of  pulses.  The  processor  will  have  the 
inherent  capability  of  operating  with  signals  of  fixed  polarization,  since  in  this  case 
all  pulse  pairs  will  be  accepted. 

The  formulations  above  are  based  on  the  assumption  that  there  are  no  echoes 
detected  in  the  coherent  receiver  which  are  aliased  in  range.  In  practice,  since 
the  Doppler  carrier  frequency  may  be  pulsed  at  approximately  1300  Hz  with  a 
corresponding  unambiguous  range  of  115  km,  this  assumption  will  not  generally  be 
valid  for  low  elevation  angles.  The  situation  is  further  complicated  in  that  these 
range-aliased  echoes  may  have  polarization  opposite  to  that  of  the  radar  antenna 
and  thus  be  detected  at  a  much  lower  power  level,  as  determined  by  the  polariza¬ 
tion  isolation  of  the  antenna.  Procedures  for  circumventing  this  problem  are  dis¬ 
cussed  in  Section  4. 


Table  2.  Percentage  of  Pulse  Pairs  Usable  for  Autocovariance 
Calculation 


(Switch  interval)  ~ 

(Doppler  Pulse  Rep.  Time) 

Autocovariance 

R(0) 

R(T  ) 

R(2t  ) 

R(3t  ) 

1 

100% 

0% 

100% 

0% 

2 

100% 

50% 

0% 

50% 

3 

100% 

67% 

33% 

0% 

4 

100% 

75% 

50% 

25% 

2.2  Absolute  and  Differential  Reflectivity 

Reflectivity  is  customarily  derived  by  averaging  the  sampled  output  of  a  log¬ 
arithmic  receiver,  in  order  to  achieve  the  widest  possible  dynamic  range.  The 
resulting  average  is  approximately  2.  5  dB  lower  than  would  be  obtained  if  the  power, 
that  is,  the  squared  signal  amplitude,  were  averaged.  For  operations  with  fixed 
polarization,  this  procedure  is  adequate  to  yield  reflectivity  estimates  with  90% 
confidence  limits  within  t  1  dB  for  typical  averaging  times. 
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A  study  of  the  statistics  of  differential  reflectivity  estimators  by  Bringi  et  al 
showed  that  the  estimates  were  consistently  more  accurate  if  the  power,  rather 
than  its  logarithm,  is  averaged.  Examples  of  their  computations  are  shown  in 
Figure  1.  The  accuracy  of  the  estimates  of  differential  reflectivity  depends  on  both 
the  number  of  samples  and  the  cross-correlation  between  signal  samples  of  opposite 
polarizations.  If  successive  pulses  are  transmitted  with  orthogonal  polarizations 
then  the  cross -correlation  is  approximately  equal  to  the  autocorrelation  of  the  re¬ 
ceived  signals  at  one -pulse  lag.  The  cross -correlation  between  successive  signals 
of  horizontal  and  vertical  polarization  is  defined  by 


phv(t)  5 


L  K(  I  AH  !  ‘ 


IE<AhAv>! 
T" 


E(  I  A, 


V  ' 


r^rr 


(29) 


where  A j j  and  A^.  are  amplitudes  of  successive  received  signals  of  horizontal  and 
vertical  polarization.  In  terms  of  autocovariance  functions  it  can  be  written  as 


phv<t)  - 


1rhv(t)! 


Irhh(0)hvv(0,I 


T7T 


(30) 


where  the  autocovariance  function  H^fT)  implies  the  summation  of  products  of 
signals  of  opposite  polarization  for  which  the  horizontally  polarized  signal  is  first, 
and  and  H^.y(0)  are  the  autocovariance  functions  computed  for  signals  of 

horizontal  and  vertical  polarizations,  respectively.  One  must  use  the  partial  sum 
of  pulse  pairs  implied  in  ^ j j \- < T >  to  avoid  the  bias  in  expected  value  resulting  from 
differential  phase  shifts  in  the  radar  and  in  the  atmospheric  propagation  medium. 

The  differential  reflectivity  is  obtained  most  accurately  from  successive  signal 
samples  with  opposite  polarization.  The  receiver  outputs  will  be  converted  to  actual 
(lower,  corrected  for  noise,  and  accumulated  in  two  registers  corresponding  to  the 
two  transmitted  polarizations.  At  the  end  of  each  averaging  interval,  the  contents 
of  one  register  are  used  to  determine  the  absolute  reflectivity  (at  horizontal  polari¬ 
zation),  and  the  estimate  of  differential  reflectivity  in  decibels  is  computed  from 
the  logarithm  of  the  ratio  of  the  two  sums.  The  cross -correlation  of  the  signal 
samples  is  to  be  estimated  from  the  Doppler  velocity  variance,  computed  from  the 
coherent  receiver  output. 


8.  Bringi,  V.N.,  Seliga,  T.  A.  ,  and  Cherry,  S.  M.  ( 1983)  Statistical  properties  of 
the  dual -polarization  differential  reflectivity  (Zjjjj)  radar  signal. 

IEEE  Trans.  Geoscience  and  Remote  Sensing,  GE-21;215-220. 
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Figure  la.  Standard  deviation  of 
estimate  of  differential  reflectivity, 
expressed  in  decibels,  as  a 
function  of  the  number  of  samples 
of  each  polarization,  based  on 
averages  of  the  squared  amplitude 
measured  at  each  polarization  for 
several  values  of  the  cross-correlation 
coefficient  between  successive  samples 
of  opposite  polarization.  [After 
Bringi  et  al  (1983),  Reference  8  4  IEEE) 


Figure  lb.  Standard  deviation  of 
estimate  of  differential  reflectivity, 
expressed  in  decibels,  as  a 
function  of  the  number  of  samples 
of  each  polarization,  for  a  given 
value  of  the  cross -correlation 
coefficient,  comparing  the  results 
due  to  averaging  power,  that  is, 
squared  amplitude,  and  the  logarithm 
of  power.  [After  Bringi  et  al  (1983), 
Reference  8  * IEEE] 


3.  SIGNAL  PROCESSOR  DESCRIPTION 

In  this  section  we  describe  the  signal  processing  hardware  and  the  computation 
of  the  algorithms  discussed  in  the  previous  section.  Current  technology  to  be  in¬ 
corporated  into  the  new  processor  makes  it  possible  to  reduce  the  overall  size  of 
the  processor  and  increase  its  computational  abilities,  in  comparison  to  the  present 
data  processor. 

Figure  2  is  a  functional  block  diagram  of  the  complete  system  showing  the 
velocity  and  reflectivity  channels.  The  velocity  channel  consists  of  a  quadrature 


receiver  with  automatic  gain  control  (ACC),  two  8 -bit  A/D  converters,  two  ground 
clutter  filters,  two  integrators  for  averaging  the  pulse  pairs  from  one  and  two-pulse 
time  lags,  two  mean  velocity  calculators,  and  a  variance  calculator  using  the 
Srivastava  algorithm  as  described  in  Section  2. 


Figure  2.  Block  diagram  of  the  radar  data  processor  illustrating  the 
calculations  of  R(t),  R(2t),  2,  and  2^ 


The  reflectivity  channel  employs  a  logarithmic  receiver,  signal  conditioner 
and  10-bit  A/D  converter,  two  separate  integrator  channels  for  averaging  the  re¬ 
turns  from  the  two  polarizations,  namely,  horizontal  (H)  and  vertical  (V),  and  a 

differential  reflectivity  calculator  (Z  ).  The  horizontal  channel  is  to  be  used  for 

L)  r\ 

absolute  reflectivity  measurements.  We  added  the  range  correction  and  range 
integrator  to  reduce  the  number  of  range  cells  for  archiving  and  display  purposes. 
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3. 1  Au  toco  variance  Calculator 


Figure  3  is  a  more  detailed  diagram  of  the  autocovariance  calculator.  For 
simplicity  only  the  single -lag  channel  is  shown  since  the  two -lag  channel  is  iden¬ 
tical.  After  the  signals  are  filtered  to  remove  ground  clutter  the  inphase  and 
quadrature  components  are  multiplied  and  summed  to  form  the  real  and  imaginary 
parts  of  the  autocovariance  given  by 


Rn=XnX(n+l)  +  Yn  Y(n+ 1) 


!n  =  Xn  Y(n+  1) 


X(n  +  D  Yn 


where  X  is  the  inphase  component,  Y  the  quadrature,  the  subscript  n  denotes  the 
pulse  delayed  by  the  interval  T,  and  n  +  1  denotes  the  undelayed  pulse.  The  two 
numbers  are  then  averaged  by  dividing  by  the  number  of  pulses  integrated  (N)  and 
summing  them  N  times  in  an  accumulator  as  shown  in  the  following  equations. 

N  R 

Re(R(T>]  =  Y,  TT  (33) 

n  =  1 


Im[R<7>]  -  £  V 
n  -  1 


The  accumulator  is  a  2048  X  25  bit  random  access  memory  (RAM).  After  accumu¬ 
lating  over  the  integration  time  the  memories  are  then  compared  range  cell  by 
range  cell  to  determine  the  larger  of  the  two  words.  The  larger  word  is  then  shifted 
in  such  a  manner  as  to  select  only  the  eight  most  significant  bits.  The  smaller 
number  is  shifted  an  identical  amount  to  keep  their  ratio  the  same.  The  eight  bits 
orovide  an  angular  resolution  of  0.  IS0  in  the  argument  of  the  autocovariance.  The 
corresponding  velocity  resolution  is  given  by 

Av  =  (0.  18°  / 180°  )v  .  (35) 


TV 


VW 


mm  "mm 


IJ.2  Average  Velocity  Calculator 

Figure  -4  is  a  simplified  block  diagram  of  the  velocity  calculator  which  com¬ 
putes  the  argument  of  the  autocovariance  on  a  range  cell  to  range  cell  basis.  The 
real  and  imaginary  numbers  from  the  accumulator  are  stored  in  the  holding  loops 
and  are  updated  each  dump  cycle.  These  holding  loops  allow  for  a  continuous  flow 
of  data  words  during  the  accumulating  time.  In  order  to  perform  the  calculations 
in  the  time  alloted  we  multiply  the  larger  number  by  the  reciprocal  of  the  smaller 
and  take  the  resultant  product  to  a  read  only  memory  (ROM)  look  up  table  to  obtain 
the  arc  tangent.  From  Figure  4  it  can  be  seen  that  any  argument  from  0  to  180'*  can 
be  calculated  unambiguously. 


(b)  45* <  9  S  90*  «'  I35*<9  1 100* 


Figure  4.  Diagram  of  mean  velocity  calculator,  illustrating  the  calculation  of 
the  argument  angle  9  in  the  domain  from  0  to  n.  The  angle  9  obtained  from 
the  arc  tangent  ROM  is  between  0“  and  45°  and  is  adjusted  to  yield  6  accord¬ 
ing  to  the  signs  and  relative  magnitudes  of  the  real  and  imaginary  parts  of 
the  autocovariance,  denoted  by  R  and  I.  Figures  a,  b,  c,  and  d  illustrate 
the  computations  when  I  >  0;  0  is  the  ones-complement  of  0  and  is  equal  to 
45°  -  0 


3.3  Doppler  Variance  Calculator 


The  Srivastava  algorithm  for  computing  the  variance  of  the  Doppler  spectrum, 

given  by  Eq,  (5),  requires  the  absolute  value  of  the  autocovariance  functions  for 

both  the  first  and  second  lag  pulse  pairs,  namely  R(t)  and  R(2r).  Since  the  aut 

covariance  function  for  each  time  lag  consists  of  two  orthogonal  vectors,  that  is, 

the  real  and  imaginary  parts,  we  take  the  square  root  of  the  sum  of  their  squares 

to  compute  the  magnitude  of  each  autocovariance  function.  Figure  5  is  a  simplified 

block  diagram  of  the  variance  calculator.  We  perform  the  high  speed  division  by 

multiplying  the  reciprocal  of  R(2t)  as  obtained  from  a  ROM.  The  natural  logarithm 

of  the  resulting  quotient,  also  derived  from  a  ROM,  is  sent  to  the  display  and 

archiving  system.  To  retain  the  greatest  versatility  in  the  processor  we  compute 
a2 

the  scaling  factor,  - ’r~’T>  by  entering  the  wavelength  (r)  and  the  pulse  repetition 

24)7T 

time  (t)  into  the  computer  that  controls  the  displays  and  performs  other  data 
processing  functions. 


Figure  5.  Details  of  the  Doppler  variance  calculator 
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3.4  Reflectivity  Calculator 


Calculation  of  differential  reflectivity,  as  discussed  in  Section  2,  requires  the 
measurement  of  the  reflectivity  factor  (Z)  at  both  horizontal  and  vertical  polariza¬ 
tions.  The  polarization  is  switched  on  a  pulse -to -pulse  basis  and  the  received 
returns  are  amplified  logarithmically  to  achieve  the  dynamic  range  needed  for 
meteorological  signals.  The  output  of  the  receiver,  which  is  the  logarithm  of  the 
instantaneous  power,  is  then  converted  to  a  10  bit  binary  number.  A  one  bit  change 
is  approximately  equal  to  a  0.  1  dB  change  in  power. 

O 

Following  the  recommendation  of  Bringi  et  al  we  convert  the  logarithmic 
output  to  power  by  means  of  a  look-up  table  in  the  form  of  a  ROM.  Floating  point 
notation  is  used  because  the  range  of  numbers  is  far  too  great  to  handle  by  straight 
binary  representation. 

Figure  6  is  a  block  diagram  of  the  reflectivity  processor.  After  the  conversion 
to  power  the  returns  are  separated  into  two  channels,  one  for  horizontal  and  the 
other  for  vertical  polarization.  An  average  is  calculated  in  each  channel  for  a 
number  of  pulses  which  is  selectable  from  2  to  1024  in  increments  of  powers  of  2. 

The  averaged  returns  are  then  stored  in  holding  loops  to  provide  a  continuous  output 
during  the  averaging  time.  After  the  holding  loops  the  average  power  in  each  channel 
is  converted  back  to  its  logarithm  by  a  ROM.  The  vertical  channel  contains  a  cor¬ 
rection  factor  for  calibration  differences  between  the  two  polarizations,  and  the 
differential  reflectivity  in  decibel  notation  is  calculated  by  taking  the  difference  of 
the  two  channels. 

The  horizontal  channel  is  chosen  for  computing  absolute  reflectivity.  We  cor¬ 
rect  low  signal -to-noise  ratios  by  subtracting  from  the  signal  an  appropriate  number 
from  a  look-up  table  which  is  based  on  the  ratio  of  the  average  signal  to  average 
noise.  This  correction  allows  for  reflectivity  values  to  be  calculated  down  to  numbers 
one  bit  above  average  noise.  After  noise  correction,  the  numbers  are  corrected  for 
range  by  adding  a  correction  which  is  based  on  the  range  cell  number  and  obtained 
from  a  ROM.  Adjacent  range  cells  can  be  averaged  to  reduce  the  number  of  range 
samples  for  archiving  purposes  and  to  provide  more  independent  samples  in  the 
averages.  The  output  is  sent  to  the  archiving  and  display  systems  for  storage  and 
real  time  display. 
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15.3  Timing  Generator 


Figure  7  shows  a  simplified  block  diagram  of  the  timing  generator.  The 
reflectivity  channel  timing  generator  and  the  velocity  dump  generator  logic  are 
shown  for  simplicity.  Upon  receipt  of  a  reflectivity  trigger  either  a  2  MHz  clock 
from  the  ground  clutter  canceller  (ACX’j )  or  the  on-board  delay  line  clock  generator 
is  divided  to  produce  either  a  1  or  2  gtsec  pulse  train.  The  delay  line  clock  genera¬ 
tor  is  used  when  the  canceller  is  not  in  use.  This  pulse  train  is  then  separated  and 
delayed  to  produce  the  various  timing  clocks  needed  for  memory  and  processor 
operation.  A  modulo  12  memory  address  counter  allows  for  the  storage  of  the 
3066  range  cells  needed  for  the  460-km  reflectivity  coverage.  With  the  use  of 
random  access  memory  in  the  processor  the  range  cells  now  span  the  entire  pulse 
repetition  interval  regardless  of  its  value. 

In  order  to  eliminate  pulse  pair  products  involving  oppositely  polarized  pulses, 
the  velocity  control  logic  was  devised.  Figure  8  shows  the  timing  diagram  for  this 
logic  for  both  the  one  and  two  lag  case  with  a  4  to  1  ratio  of  velocity  to  reflectivity 
PRF.  The  shaded  area  denotes  the  time  for  processing.  At  all  other  times  the 
output  of  the  complex  multiplier  is  set  to  zero  and  its  output  is  not  accumulated. 

Only  the  processed  triggers  are  counted  to  determine  the  integration  time  of  the 
processor.  The  number  of  pulses  integrated  is  accomplished  by  selecting  the 
proper  output  for  the  counter.  It  is  selectable  from  2  to  1024  in  intervals  of  powers 
of  2. 


4.  OPERATIONAL  PROCEDURES  AND  PARAMETERS 

Outputs  of  the  data  processor  will  be  displayed  and  archived  by  a  computer 
controlled  data  system.  We  anticipate  that  the  absolute  reflectivity,  differen¬ 
tial  reflectivity,  floppier  mean  velocity,  and  Doppler  variance  will  be  routinely 
available  for  display  in  plan  position  or  range-height  coordinates.  Derived 
parameters  such  as  number  concentration  and  drop  size,  rainfall  rate,  wind 
gradients,  and  turbulence  intensity  will  be  calculated  by  the  computer. 

Opera'ional  parameters  for  the  radar  will  include  the  pulse  repetition  frequency 
(PRF)  of  each  of  the  two  transmitters  and  the  polarization  switching  frequency.  It 
is  expected  that  the  polarization  will  be  switched  at  the  PRF  of  the  transmitted 
signal  used  for  reflectivity  measurement  and  that  the  PRF  of  the  transmitted  signal 
used  for  Doppler  velocity  measurement  will  be  3  or  4  times  higher,  as  discussed  in 
Section  2.  Switch -selectable  parameters  for  the  data  processor  include  the  number 
of  pulses  to  be  averaged  in  time  and  in  range  for  the  computation  of  reflectivity,  the 
number  of  pulses  to  be  averaged  for  input  to  the  automatic  gain  control,  and  the 
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number  of  pulse  pairs  to  be  averaged  for  the  autoeovariance  estimates.  The  radar 
wavelength  and  PRF  of  the  Doppler  carrier  frequency  will  be  input  through  the  com¬ 
puter  for  use  in  scaling  the  velocity  and  variance  displays,  and  the  peak  power  of 
the  transmitted  signal  used  for  reflectivity  measurement  will  be  input  for  the  calcu¬ 
lation  of  reflectivity.  All  these  parameters  will  be  archived,  along  with  the  primary 
data  processor  outputs  and  ancillary  information  such  as  date,  time,  and  location. 
The  ranges  and  increments  of  the  control  parameters  are  shown  in  Table  3. 


H _ | _ V _ | _ H 

Reflectivity  channel  trigger 


Velocity  channel  triaaer 
One-pulte  lag  control 

1  _ . Wm  Wm  . 

Two-pulse  lag  control 

Figure  8.  Timing  diagram  for  the  computation  of 
differential  reflectivity  and  autocovariance  for  the 
case  in  which  the  Doppler  carrier  PRF  is  four 
times  the  reflectivity  carrier  PRF.  The  control 
signals  (third  and  fourth  lines)  are  used  to  suppress 
the  complex  multiplier  outputs  that  correspond  to 
pairs  of  pulses  with  opposite  polarizations 


Table  3.  Data  Processor  Characteristics 


Feature 

Specification 

Range  Cells 

Number 

Dependent  on  the  pulse  repetition  time. 

Spacing 

1  or  2  fisec,  selectable 

Doppler  Channel: 
Integration  Type 

Pulses  Integrated 

Rectangular 

2n,  n  selectable  1  to  10 

Range  Cells  Integrated 

Moving  average  of  2n,  n  selectable  0  to  3 

Reflectivity  Channel: 
Integration  Type 

Pulses  Integrated 

Rectangular 

2n,  n  selectable  1  to  10 

Range  Cells  Integrated 

Moving  average  of  2n,  n  selectable  0  to  3 

AGC  Channel: 

Integration  Type 

Exponential 

Pulses  Integrated 

2n,  n  selectable  0  to  10 

Recording: 

Magnetic  Tape 

Analog,  14  track,  serial  format 

Video  Display: 

Type 

Real-time,  programmable,  color 
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Several  procedures  may  be  used  to  circumvent  the  problem  of  range- aliased 
echoes  which  may  be  present  in  the  Doppler  channel.  Each  of  these  involves  a 
trade-off  of  sampling  time,  time  resolution  of  the  averaged  data,  or  numerical 
accuracy.  Lf  the  PRF  of  the  Doppler  carrier  frequency  is  decreased  then  the 
increased  unambiguous  range  will  encompass  all  echoes,  but  the  maximum  unambig¬ 
uous  Doppler  velocity  will  be  decreased.  Furthermore,  the  increased  time  interval 
between  pluses  of  the  reflectivity  carrier  frequency  will  yield  a  decreased  cross¬ 
correlation  between  successive  samples  >f  reflectivity  at  orthogonal  polarizations, 
leading  to  reduced  accuracy  of  the  differential  reflectivity.  A  second  procedure 
is  to  operate  sequentially  with  fixed  and  switched  polarization,  obtaining  reflectivity 
and  Doppler  velocity  measurements  from  the  signals  of  fixed  polarization  and  ob¬ 
taining  absolute  and  differential  reflectivity  from  the  signals  transmitted  with  al¬ 
ternating  polarization.  This  procedure  would  necessitate  double  scanning  at  low- 
elevation  angles  and  thus  increase  the  time  between  successive  observations  of  any 
given  region  of  space.  A  third  alternative  is  to  transmit  more  than  four  pulses  of 
the  Doppler  carrier  frequency  at  each  polarization  setting.  If  the  returned  signal 
is  sampled  beginning  with  the  third  pulse  repetition  interval,  then  aliased  echoes 
of  the  unwanted  polarization  will  not  be  present.  Although  this  procedure  permits 
the  high  unambiguous  velocity  interval  to  be  maintained,  the  numerical  accuracy 
of  the  differential  reflectivity  will  be  degraded  because  of  the  longer  time  between 
the  sampling  of  orthogonal  polarizations. 


5.  CONCLUSIONS 

We  have  developed  a  design  for  a  -lata  processor  capable  of  computing  absolute 
reflectivity,  differential  reflectivity,  Doppler  mean  velocity,  and  Doppler  variance 
from  received  signals  in  a  radar  system  transmitting  pulses  of  alternating  horizontal 
and  vertical  polarization.  The  velocity  parameters  are  to  be  derived  from  estimates 
of  the  signal  autocovariance  function  computed  by  summing  the  complex  products  of 
signal  amplitudes  at  zero,  one,  and  two-pulse  time  lags.  The  key  design  factor  here 
is  the  need  to  delete  from  the  summation  the  products  of  signals  transmitted 
with  different  polarization.  This  need  arises  from  the  effects  of  differential  phase 
shift  between  the  two  polarizations  in  propagation.  The  reflectivity  processor  uses 
floating-point  computational  hardware  to  average  the  received  power,  rather  than  its 
logarithm.  This  design  feature  permits  the  most  accurate  computation  of  the  differ¬ 
ential  reflectivity.  The  processor  will  be  used  with  a  computer-controlled  display 
and  archiving  system  to  produce  real-time  displays  of  the  primary  processor 
outputs  and  derived  meteorological  quantities. 
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